Baryon Isocurvature Perturbation in the Affleck-Dine Baryogenesis mechanism 
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We propose a primeval baryon isocurvature model in the context of the Affleck-Dine (AD) 
baryogenesis mechanism. Quantum fluctuations of the AD field during inflation produce fluctuations 
in the baryon number. We obtain a model which gives both the appropriate baryon to entropy ratio 
of the order lO" 11 and the gaussian baryon isocurvature perturbation of the order 10 -3 at the break 
of the spectrum IMpc with the steep spectrum n = —1.5 on large scales > IMpc. 
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The primeval baryon isocurvature model (PBI) which 
attempts to account for the large scale structure of the 
universe only by the baryon and the radiation has sev- 
eral appealing properties which are lacking in standard 
CDM model [fij. Under the assumption that there ex- 
ist primeval fluctuations in the baryon to entropy ratio 
with a steep spectrum on large scales, PBI fits two sets 
of observations ||[|. First, large power on small scales in 
the fluctuations spectrum produces galaxy at a relatively 
early stage, which is supported by recent galaxy surveys. 
Second, the large scale correlations in the distributions 
of galaxies would be produced at matter-radiation Jeans 
scale ~ lOOMpc. 

So far, most of studies about PBI are started by as- 
suming the existence of ad hoc primordial fluctuations. 
To realize PBI in the physics of early universe, one needs 
a model of baryogenesis which can give both the appro- 
priate baryon to entropy ratio ub/s of the order 10 
and its fluctuations with the steep spectrum. The tilted 
spectrum is necessary to suppress the cosmic microwave 
background (CMB) anisotropy still keeping sufficient am- 
plitude on the galaxy scale A possible model is the 
one based on the grand unified theory (GUT) baryogen- 
esis and inflation. In Ref ||, the baryogenesis from the 
decay of the heavy Majorana lepton with both hard and 
soft CP violation is considered and the non-gaussian fluc- 
tuations with the steep spectrum are obtained. Gener- 
ally, however, the GUT scale (Mq) inflation produces the 
adiabatic perturbation of the order (Mc/M p ) 2 ~ 10~ 5 , 
where M p is the planck scale, so one needs the artificial 
suppression to ensure it's subdominance to the isocurva- 
ture perturbation. If one takes the lower energy as the 
scale of the inflation, the adiabatic perturbation is sup- 
pressed, but the reheating temperature also decreases, so 
the GUT baryogenesis would not work well. Then, we 
need the model of baryogenesis at lower energy scales. 
Unfortunately, the weak scale baryogenesis in the stan- 
dard model (SM) of the particle physics is unable to 
produce the sufficient baryon asymmetry ||, but a vi- 
able alternative exists in supersymmetric extension of 
SM, i.e., the Affleck-Dine (AD) mechanism The AD 
mechanism produces the baryon asymmetry efficiently 
and it is known that scale invariant fluctuations in the 



baryon number are induced by the quantum fluctuations 
of the field responsible for the baryogenesis (AD field) 
during inflation 0]. However, as mentioned above, the 
steep spectrum is crucial to realize PBI. In this letter, 
we will show the tilted spectrum of fluctuations is imple- 
mented in AD baryogenesis mechanism by coupling the 
AD field to another field Q. 

In the AD mechanism, the baryogenesis occurs from 
the coherent condensation of the AD field, which carries 
B — L U(l) charge, where B and L are the baryon num- 
ber and the lepton number. The potential of the AD field 
is given by [§ 
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where TO3/2 is the gravitino mass, H is the Hubble pa- 
rameter, M is some large mass scale, A is the dimension- 
less coupling constant, n is some integer larger than 4, 
a — ae l9a and A — Ae l8A are some complex numbers, and 
c is a real constant. Here we will consider the case n = 4 
and A = 1 Jll| . The important point of the model of Ref 
S is that the negative effective mass term is assumed 
to exist during inflation. During inflation, the AD field 
has large value and subsequent to inflation, it spirals into 
the origin exciting the baryon number and decays to the 
baryon after reheating. Following, we will describe how 
baryon isocurvature fluctuations with steep spectrum can 
be obtained in this model of bayogenesis. 

(1) During inflation (H 3> W3/2)) we can ignore the 
terms proportional to m 3 / 2 . Due to the negative mass 
term proportional to H, the potential has global mini- 
mum far from the origin and undergoes quantum fluctu- 
ations of the order H about the minimum. || . We con- 
sider the scales which are still out of horizon when the 
AD field decays and creates baryons, so the spectrum 
shape of baryon number fluctuations is determined by 
that of the AD field during inflation. We want the tilted 
spectrum favor to small scales in order to suppress the 
CMB anisotropy [Q . Generally, the tilt of the spectrum 
can be obtained if the modes of the fluctuations have 
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the mass proportional to the Hubble parameter. This is 
because amplitudes of modes which are massive decrease 
with time even after they cross the horizon. For example, 
if the mode has the mass m 2 = m 2 H 2 and the Hubble 
parameter Hi during inflation is constant, the mode am- 
plitudes decrease as R~ 3 / 2 + y (y = (9/4 — ml) 1 / 2 ) after 
the mode crosses the horizon. The mode amplitudes are 
Hi when it crosses the horizon, so at the end of the in- 
flation, the mode function becomes 
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where Rk and R e are scale factors when the mode ex- 
its the horizon and the inflation ends. The tilt from the 
scale invariant spectrum n — — 3 is given by An = 3 — 2v, 
which is larger as the mass becomes large. This could 
produce the isocurvature perturbation with the tilted 
spectrum fl2| . However, since the comoving scale of 
galaxies is much larger than the horizon scale at the 
inflationary epoch (k ga i axy /R e <C Hi), the fluctuation 
amplitude on the galaxy scale suffers from severe sup- 
pression To avoid this problem we must change the 
break of the spectrum from Hi to another scale. This 
can be done if (p couples to the other field ip which gives 
the mass m 2 oc ip 2 to (p. First, we assume the power-low 
inflation in order to obtain the desired behavior of ip [^0) . 
The potential of the inflaton \ is 



U(\) = A 4 exp 




(3) 



where A is some mass scale and e, which is less than unity, 
determines the expansion rate during inflation. The scale 
factor and the Hubble parameter are R oc t l l\ H = 
l/(et). Second, we couple the AD field to the other real 
field ip and assume the negative coupling. The potential 
is 

V(<P, = VM) + \^ + i/iV - + 0* 2 ), (4) 

where 7 and (3 are some positive coupling constants which 
we assume tiny (~ m 3 / 2 /M) and (3ip 2 , fi 2 ^> "f\(p\ 2 Q. In 
this case, we can ignore the effect of coupling to <p in the 
evolution of ip. Initially ip is assumed to have the large 
value then at early times the self-coupling term is dom- 
inant, and ip slowly decreases with t (oc t^ 1 ) so evolves 
like H (ip= (e/^ 1 / 2 )(3e- 1 - 2f/ 2 H). The behavior of iP 
changes at t\ when the mass term and the self-coupling 
term become comparable. At that time ip = /z/3 -1 / 2 , so 
t\ is determined by the mass [i (t\ = [sT 1 ^ 1 / 2 ). Later, 
the mass term becomes dominant and the amplitude of 
ip decays rapidly (ip oc i~ 3 / 2e ) |[o]]. Following, for sim- 
plicity, we assume a = 0. For t < t±, due to the coupling 
with ip, the minimum of the AD field <p — <pe l6 and the 
mass around the minimum m 2 = m 2 H 2 is given by 

Ooeff =0,m 2 ee// = 4 £ (3-2e)( 7 //3), 
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where c e // = c+ 2e(3 — 2e)(p/ / (3) and m^e// is the mass 
of the amplitude and mg e ff is that of the phase. For 
t > ii, ip decays, so the minimum and the mass are given 
by 
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For t < t\, the curvature around the minimum is large 
so fluctuations are massive (w e // > m). Then, the am- 
plitudes of the modes which exit the horizon for t < t\ 
decrease with time more rapidly after crossing the hori- 
zon than the modes which exit the horizon for t > t±. So, 
the decay of ip for t > t\ gives the break to the spectrum 
at k\ which is the wavenumber of the mode which exits 
the horizon at t = t\. The modes of k < k\ have the 
steeper spectrum than the modes of k > k\ . Taking into 
consideration that the Hubble parameter varies with time 
Jbij , we obtain the mode functions of the modes k > k\ 
at the end of the inflation 
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where 
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and Hi and Ri are the Hubble parameter and the scale 
factor when the mode with wavenumber k\ exits the hori- 
zon. The mode functions of the modes k < k\ are ob- 
tained by replacing vg to ug e f f which is given by replacing 
m*g to m*g e ff Notice that the break of the spectrum is 
ki rather than the Hubble parameter during inflation as 
required. The spectrum index n = —2vg e ff for k < k\ 
is larger than n = —2vg for k > k\. Fluctuations of the 
amplitude are given by replacing qg, vg to v$ respec- 
tively. 

(2)Subsequent to inflation (H > to 3 / 2 ), we assume the 
universe enters the matter era during inflaton oscilla- 
tion, so the Hubble parameter varies with time. For 
H > TO3/2, we can still ignore the terms proportional 
to m 3 / 2 - Since we assume a = 0, the phase and its fluc- 
tuations are constant (8 = 9q, 59 = 59\h,,)- In the case 
c is not so tiny the field tracks the fixed point near the 
instantaneous minimum of the potential and decreases 
with time (~ (HM) 1 / 2 oc t" 1 / 2 ) [§ and the fluctuations 
vary in the same way (oc i -1 / 2 ), thus during this times 
8(p/(po = const.. If c is very tiny or vanishes Jig] , the 
higher terms becomes important, so the AD fields is near 
critically damped (V"(<p>) ~ H 2 ). In this case there also 
exists the attracting point (~ (HM) 1 / 2 ), so the situation 
is almost the same with the former case [HI. 
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(3) At H ~ W3/2, the low energy supersymmetry break- 
ing terms proportional to m 3 / 2 become important and 
the B — L non-conserving A-term arises. The AD field 
falls toward the origin. The evolution of the AD field at 
this epoch can be considered as the motion of a point- like 
particle in the two dimensional plane (4>r, </>/), where 4>r 
and 4>i are the real and imaginary part of the AD field. 
The baryon charge n B = 2\4>\ 2 9 = 2(cf>j(j)R — cf)R<f)i) corre- 
sponds to the angular momentum of the particle, and its 
non-conservation is caused by [/(l)-asymmetric A-term 
p6| . In the case A = c = r\ <C 1, we can obtain the 
approximate solution. The ansatz of the solution is 



= ^ur sin(i - U + (3 R ( 



1 



t U 

h = ja/ sm(t -U+ /3i(- - -!- 



) + S R ), 

)+Si), (8) 



where we rescale <f> — ► (m 3 / 2 M ) 1//2 </> and t — > m 3 ±t. The 
initial condition is set at t = tj. Until t = ti, we as- 
sume the AD fields tracks the fixed point, so at that time 
4> = 4>i, 4> — —{\/2ti)4>, 9 — 9 , 9 — 0. 4>i is determined 
by the time when the baryon number is excited (t = ti). 
At this time H ~ 1713/2, so 4>i is ~ 0(1). We obtain the 
solution to the first order of 77 
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and £r = <5j = arctan(2), where </>fli = <f>i cos(#a/4), 4>h = 
0jSin(#A/4). Here we neglect the terms proportional to 
t p (p < —2) and the high frequency terms. The baryon 
number at late times tends to 

^~§(f) (mj /2 M)Afr 4 sme A . (10) 

We have done numerical calculations and find that this 
approximate solution agrees with numerical results up 
to a factor 2. We consider fluctuations of the field as 
fluctuations in the initial condition of this motion. For 
H > TO3/2J 89 and 64>/4>o are approximately constant, so 
fluctuations in the baryon number become 
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H = H e 



In the case a = 0, c 7^ 0, fluctuations of the amplitude are 
suppressed by the factor (i? e /i?i) _91 * compared to those 
of the phase so negligible. For H < 7713/2 A-term becomes 
negligible, so the baryon number and its fluctuations are 
constant in the comoving volume. 

(A) After reheating, the inflaton \ converts to the ra- 
diation and the entropy is produced. The fractional 
energy in the AD condensate at this time is p<f,/pi ~ 



(m 3 / 2 M)/Mp <C 1, so the entropy is dominated by in- 
flaton decay. To ensure the condensate is not damped 
by thermal scattering before the baryon asymmetry is es- 
tablished, the reheating temperature should be Tr < 10 6 
Gev S. For the reheating temperature above the weak 
scale, B — L is converted to the baryon number. 

We will show that both the appropriate baryon to en- 
tropy ratio and its fluctuations with the appropriate am- 
plitude and spectrum shape are obtained by choosing 
specific parameters. The baryon to entropy ratio is fl 



n B 
s 



75 
G4 



A$sin0 A lO- 10 ' Tr 



Vl0 6 Gev 



10- 3 M 
M„ 



(12) 



When we choose A = c = 0.1, & = 1,0a = 7r/4, M = 
10 3 M P , and Tr = 10 6 Gev, this becomes n B /s = 1CT 11 , 
which is consistent with the constraint from primordial 
nucleosynthesis. The power spectrum of the entropy per- 
turbation is given by 
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where 
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We set the expansion rate during inflation at e = 0.28. 
The rest parameters are determined so as to satisfy as- 
trophysical constraints. First, in order to suppress the 
CMB anisotropy, the spectrum index should be n > — 1.5 
on large scales [||J|. So, we choose 7//? = 0.57, which 
corresponds to rvQ Be ^ = 1.56, then the spectrum index 
becomes n — —2h>g e ff = —1.5. The amplitude and the 
break are determined by J3 normalization || 
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where x\ = 2ir/ki is the comoving scale of the break, 
Do/D rec is the growth factor from recombination to 
present and P rec is the power spectrum at recombina- 
tion; P rec = P B (fc>fe eg ), P rec = P B (k/k eq ) 4 (fc<fc eg ). 
Here k eq is the comoving wavenumber corresponding to 
the horizon scale at the time of equality between baryon 
and radiation density. For example, we fix the break at 
xi = IMpc H, then the amplitude Q should be 10~ 3 . 
To obtain the appropriate baryon to entropy ratio, H e 
should be larger than m 3 / 2 - So we take H e = 10 5 Gev 
and Hi = 10 12 Gev, then Q ~ 10 -3 . The expansion fac- 
tor during inflation is given by 
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This is appropriate to adjust the scale of k\ to IMpc. 
The value of the Hubble parameter during inflation is 
relatively low, so the primeval adiabatic perturbation is 
naturally subdominant. In fact, the curvature perturba- 
tion 1Z of the inflaton \ can be estimated by the standard 
formula; 
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which is negligible compared to the baryon isocurvature 
perturbation. We make assumption that during inflation 
the coupling to the AD field does not affect the evolution 
of tp (/3ip 2 ,[i 2 ^> 7|^| 2 )- Since the mass of ip is related 
to H x by fi 2 = eHl and i\) 2 > eft" 1 !! 2 for t < h, this 
condition is satisfied if 7 < eHi/M, which is consintent 
with the assumption 7 ~ j3 ~ 

In summary, we have obtained the PBI, which gives 
both the appropriate baryon to entropy ratio ~ 10~ n 
and the gaussian baryon number perturbation of the or- 
der ~ 10~ 3 at the break x\ — IMpc with the steep spec- 
trum n = —1.5 on large scales x > x\. The inflation 
occurs at relatively low scale such as 10 12 Gev, so the 
adiabatic perturbation generated by the inflaton is neg- 
ligible naturally. The e-folding number during inflation 
is sufficient to solve the horizon problem. The reheat- 
ing temperature is higher than the weak scale and low 
enough to avoid the gravitino problem. 

The most essential assumption of our model is the cou- 
pling between AD field and the field ip. The tiny coupling 
constant 7 ~ m 3 / 2 /M is needed |l3|. In our model, the 
coupling to tp produces the large initial values for baryo- 
genesis in early stage of inflation, so even in the case 
when the terms coming from finite energy supersymme- 
try breaking are small or vanishes [p^|, our model still 
works well. The exact value of coefficients of the poten- 
tial relevant to the AD field and vb may be determined 
by high energy physics one has not known well and a set 
of parameters we take here is not unique. Our purpose 
in this letter is to show there exists at least one PBI in 
physical model of baryogenesis. 
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